Measurement of cross sections of exclusive e + e — > VP processes 

at = 10.58 GeV 

(Belle Collaboration) 

K. Belous 12 , M. Shapkin 12 , I. Adachi 9 , H. Aihara 43 , K. Arinstein 1 ' 32 , V. Aulchenko 1 * 32 , A. M. Bakich 39 , 
V. Balagura 14 , E. Barberio 22 , W. Bartel 9 , A. Bay 19 , M. Bischofberger 24 , A. Bondar 1 * 32 , A. Bozek 28 , 
M. Bracko 21 ' 15 , T. E. Browder 8 , P. Chang 27 , Y. Chao 27 , A. Chen 25 , B. G. Cheon 7 , I.-S. Cho 47 , 
S.-K. Choi 6 , Y. Choi 38 , J. Dalseno 9 , M. Dash 46 , A. Drutskoy 3 , S. Eidelman 1 - 32 , D. Epifanov 1 ' 32 , 
N. Gabyshev 1 ' 32 , A. Garmash 1 - 32 , H. Ha 17 , Y. Horii 42 , Y. Hoshi 41 , W.-S. Hou 27 , H. J. Hyun 18 , T. Iijima 23 , 
K. Inami 23 , A. Ishikawa 35 , R. Itoh 9 , M. Iwasaki 43 , N. J. Joshi 40 , D. H. Kah 18 , N. Katayama 9 , H. Kawai 2 , 
T. Kawasaki 30 , H. O. Kim 18 , J. H. Kim 38 , Y. I. Kim 18 , Y. J. Kim 5 , B. R. Ko 17 , P. Krizan 20 ' 15 , 
P. Krokovny 9 , R. Kumar 34 , A. Kuzmin 1 - 32 , Y.-J. Kwon 47 , S.-H. Kyeong 47 , M. J. Lee 37 , S.-H. Lee 17 , 
T. Lesiak 28 ' 4 , A. Limosani 22 , C. Liu 36 , D. Liventsev 14 , R. Louvot 19 , A. Matyja 28 , S. McOnie 39 , 
H. Miyata 30 , R. Mizuk 14 , T. Mori 23 , Y. Nagasaka 10 , S. Nishida 9 , O. Nitoh 45 , T. Ohshima 23 , S. Okuno 16 , 

H. Ozaki 9 , P. Pakhlov 14 , G. Pakhlova 14 , C. W. Park 38 , H. K. Park 18 , R. Pestotnik 15 , L. E. Piilonen 46 , 
A. Poluektov 1 ' 32 , Y. Sakai 9 , O. Schneider 19 , C. Schwanda 13 , K. Senyo 23 , M. E. Sevior 22 , V. Shebalin 1 - 32 , 

C. P. Shen 8 , J.-G. Shiu 27 , B. Shwartz 1 ' 32 , A. Sokolov 12 , S. Stanic 31 , M. Staric 15 , J. Stypula 28 , 
T. Sumiyoshi 44 , G. N. Taylor 22 , Y. Teramoto 33 , K. Trabelsi 9 , S. Uehara 9 , Y. Unno 7 , S. Uno 9 , Y. Usov 1 ' 32 , 
G. Varner 8 , K. E. Varvell 39 , K. Vervink 19 , A. Vinokurova 1 - 32 , C. H. Wang 26 , P. Wang 11 , Y. Watanabe 16 , 
R. Wedd 22 , E. Won 17 , B. D. Yabsley 39 , H. Yamamoto 42 , Y. Yamashita 29 , V. Zhilich 1 - 32 , V. Zhulanov 1 * 32 , 

T. Zivko 15 , A. Zupanc 15 , O. Zyukova 1 - 32 

1 Budker Institute of Nuclear Physics, Novosibirsk, Russian Federation 
2 Chiba University, Chiba, Japan 
8 University of Cincinnati, Cincinnati, OH, USA 
4 T. Kosciuszko Cracow University of Technology, Krakow, Poland 
5 The Graduate University for Advanced Studies, Hayama, Japan 
6 Gyeongsang National University, Chinju, South Korea 
7 Hanyang University, Seoul, South Korea 
8 University of Hawaii, Honolulu, HI, USA 
9 High Energy Accelerator Research Organization (KEK), Tsukuba, Japan 
10 Hiroshima Institute of Technology, Hiroshima, Japan 
11 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing, PR China 
1 11 Institute for High Energy Physics, Protvino, Russian Federation 
18 Institute of High Energy Physics, Vienna, Austria 
1J > Institute for Theoretical and Experimental Physics, Moscow, Russian Federation 
15 J. Stefan Institute, Ljubljana, Slovenia 
16 Kanagawa University, Yokohama, Japan 
17 Korea University, Seoul, South Korea 
18 Kyungpook National University, Taegu, South Korea 
1!> Ecole Polytechnique Federate de Lausanne, EPFL, Lausanne, Switzerland 
20 Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana, Slovenia 
21 University of Maribor, Maribor, Slovenia 
22 University of Melbourne, Victoria, Australia 
28 Nagoya University, Nagoya, Japan 
2 ^Nara Women's University, Nara, Japan 
25 National Central University, Chung-li, Taiwan 
2e National United University, Miao Li, Taiwan 
27 Department of Physics, National Taiwan University, Taipei, Taiwan 
28 H. Niewodniczanski Institute of Nuclear Physics, Krakow, Poland 
'^Nippon Dental University, Niigata, Japan 
80 Niigata University, Niigata, Japan 
8 University of Nova Gorica, Nova Gorica, Slovenia 
82 Novosibirsk State University, Novosibirsk, Russian Federation 
88 Osaka City University, Osaka, Japan 



Preprint submitted to Physics Letters B 



November 9, 2009 



3J i Panjab University, Chandigarh, India 
35 Saga University, Saga, Japan 
se University of Science and Technology of China, Hefei, PR China 
s7 Seoul National University, Seoul, South Korea 
ss Sungkyunkwan University, Suwon, South Korea 
s9 University of Sydney, Sydney, NSW, Australia 
Jf0 Tata Institute of Fundamental Research, Mumbai, India 
Jfl Tohoku Gakuin University, Tagajo, Japan 
^Tohoku University, Sendai, Japan 
J < '"' Department of Physics, University of Tokyo, Tokyo, Japan 
*■* Tokyo Metropolitan University, Tokyo, Japan 
4 5 Tokyo University of Agriculture and Technology, Tokyo, Japan 
4 6 IPNAS, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA 
^ Yonsei University, Seoul, South Korea 



Abstract 

The cross sections for the reactions e + e~ — > ^77, <prj' ', prj, prj have been measured using a data sample of 
516 fb -1 collected with the Belle detector at the KEKB asymmetric-energy e + e~ collider. The corresponding 
values of the cross sections are: 1.4 ± 0.4 ± 0.1 fb (0r?), 5.3 ± 1.1 ± 0.4 fb (<jrrf), 3.1 ± 0.5 ± 0.1 fb {prj) and 
3.3 ± 0.6 ± 0.2 fb (prj'). The energy dependence of the cross sections is presented using Belle measurements 
together with those of CLEO and BaBar. 



High-statistics data samples accumulated at B- 
factories allow a study of rare exclusive two-body 
processes in e + e~ annihilation. An example of 
such a highly suppressed process is the reaction 
e + e~ — > VP, where V and P stand for Vec- 
tor and Pseudoscalar mesons, respectively. It has 
been observed [l|, 0] that double charm produc- 
tion in e + e~ — > J/ipVc has an unexpectedly high 
cross section. The basic diagram for double charm 
production is very similar to the one describing 
e + e~ — > (prj{rf) where c quarks are replaced by s 
quarks. Thus, comparison of the two reactions may 
contribute to better understanding of the underly- 
ing physics. In addition, we investigate the pro- 
cesses e + e~ — > prj(rj') 1 which also belong to the VP 
class with a different isospin configuration and light 
quarks only. 

Some of the e + e~ — > VP reactions have pre- 
viously been measured at different center-of-mass 
(CM) energies: in the DM1 experiment at y/s be- 
tween 1.4 and 2.18 GeV [3], in the CLEO experi- 
ment at y/s = 3.67 GeV [4[ and by the BaBar col- 
laboration at y/s between 1 and 3 GeV using initial- 
state radiation (ISR) [| and at y/s = 10.58 GeV Q. 
The cross section of the process e + e~ — > prj was 
also measured by the BES collaboration at y/s = 
3.65 GeV Q- The reaction e + e~ — > (jnf has not yet 
been observed, nevertheless the upper limit on its 



cross section set by CLEO [4] can be useful for dis- 
crimination between models that predict different 
energy dependences. 

The QCD-based models predict the energy de- 
pendence for the process e + e~ — > VP to be 
1/s 4 @, Q while the cross section for the process 
e + e~ — > (f>r/ measured by CLEO and BaBar favors 
a 1/s 3 dependence. The form factor for the pro- 
cess e + e~ — > VP is expected to have a 1/s de- 
pendence [Io| . Recently theoretical calculations of 
e+e~ — * VP cross sections have been published, 



which use the light cone approach 11 



Pre- 



dictions are given for two values of y/s: 3.67 and 
10.58 GeV. The authors of Ref. [ll| claim that 
their results favor a 1/s 3 dependence. In Ref. [13] 
a r~j 1/s 4 is expected in the limit s — > 00. 

The analysis presented here is based on data 
taken at the T(4S) (y/s = 10.58 GeV) with the 
Belle detector at the KEKB asymmetric-energy 
e + e~ collider [lj|. The total integrated luminosity 
of the on-resonance sample used in the analysis is 
516 fb -1 . To check that the processes e + e~ — > VP 
are due to a single-photon annihilation and that the 
hadronic decay of the T(45) does not give a signifi- 
cant contribution, we use a 58 fb -1 data sample col- 
lected 60 MeV below the resonance peak. All the 
observed off-resonance signals are consistent with 
those at the T(45) resonance within statistical cr- 



rors. We use these data to set upper limits for the 
branching ratios of T(45) — > VP. 

A detailed description of the Belle detector is 
given elsewhere 14j. We mention here only the de- 
tector components essential for the present analy- 
sis. Charged tracks are reconstructed from hit in- 
formation in the central drift chamber (CDC) lo- 
cated in a 1.5 T solenoidal magnetic field. Trajecto- 
ries of charged particles near the collision point are 
provided by a silicon vertex detector (SVD). Pho- 
ton detection and energy measurements are per- 
formed with a CsI(Tl) electromagnetic calorimeter 
(ECL). Identification of charged particles is based 
on the information from the time-of-flight coun- 
ters (TOF) and silica aerogel Cherenkov counters 
(ACC). The ACC provides good separation be- 
tween kaons and pions or muons at momenta above 
1.2 GeV/c. The TOF system consists of 128 plastic 
scintillation counters and is effective in K/tt sepa- 
ration for tracks with momenta below 1.2GeV/c. 
Low energy kaons are also identified using specific 
ionization (dE/dx) measurements in the CDC. 

In order to identify hadrons, for each of the three 
hadron types i (i = 7r, K and p) a likelihood Li 
is formed using information from the ACC, TOF, 
and dE/dx measurements from the CDC. Kaons are 
selected with the requirement Lk /(Lk + L„) > 0.6, 
which has an efficiency of 90% and 6% probability 
to misidentify a pion as a kaon. All charged tracks 
that are not identified as kaons are considered to be 
pions. 

Signal candidates are selected in two steps. Ini- 
tially, events with low multiplicity are selected by 
requiring that the number of charged tracks in an 
event be two or four with zero net charge and each 
track have a momentum transverse to the beam axis 
(pt) larger than 0.1 GeV/c; and that each track ex- 
trapolate to the interaction point (IP) within 1 cm 
transversely and within 5 cm along the beam di- 
rection. To suppress background from Bhabha and 
events, the sum of the absolute values of mo- 
menta of the first and second highest momentum 
tracks are required to be less than 9GeV/c. At 
least one track with p t above 0.5 GeV/c is required. 
Beam-related background is suppressed by requir- 
ing that the position of the reconstructed event ver- 
tex be less than 0.5 cm from the IP in the transverse 
direction and less than 3 cm from the IP along the 
beam direction. 

Photons are defined as ECL clusters with energy 
deposits above 200 MeV that are not associated 
with charged tracks. Neutral pion candidates are 



formed from pairs of photons with invariant masses 
in the range 120 to 150MeV/c 2 . The two-photon 
invariant mass resolution in the mass region of the 
7T° is about 6MeV/c 2 . For the 77—5-77 reconstruc- 
tion we use only photons that do not form a ir° 
candidate with any other photon. The invariant 
mass of 77 — ► 77 candidates should lie in the range 
0.5 — 0.6GeV/c 2 . The two-photon invariant mass 
resolution in the mass region of the 77 meson is about 
20MeV/c 2 . 

After preselection we apply the following require- 
ments to extract the exclusive VP final states: 

• The difference between the energy of VP can- 
didates in the CM frame and y/s of KEKB 
should be between -0.3 and +0.2 GeV0 

• The angle between V and P candidates in the 
CM frame should be larger than 175 degrees. 

We consider the following decay modes of vector 
and pseudoscalar mesons: 0(1020) — * K + K~ , p — > 
7r + 7r~, 77 — * 77, 77' — * 7r + 7r~7, 77' — ► r)Tr + Tr~. After 
the application of the requirements listed above we 
observe significant concentrations of events in the 
scatter plots near the masses of corresponding vec- 
tor and pseudoscalar mesons. These scatter plots 
are shown in Figs. HJa-f). 

A two-dimensional unbinned likelihood fit is ap- 
plied to extract the signal yields for the above reac- 
tions. We assume the mass distributions of vector 
and pseudoscalar particles to be uncorrelated; thus 
the distributions in the scatter plots of Fig. 1 can be 
represented as the product of two one-dimensional 
probability density functions (PDF), one for each 
dimensionH To fit the 0(1020) -> K+K~ and 
p — > 7r + 7r~ invariant mass distributions we use a 
non-relativistic Breit-Wigner function. The 77 — > 77 
and 77' — > 7r + 7r _ 7, i]Tr + ir~ invariant mass distri- 
butions are fitted with Gaussians. The background 
for the K + K~ system is described by the prod- 
uct of the threshold function (m(K + K~) — Too)" 
and a first-order polynomial, where a is a free 
parameter and too = 2to^-+ . The fit range for 
m(K+K~) extends from to to 1.12GeV/c 2 . The 



This cut allows for low energy photon radiation up to 
0.3 GeV to be present. 

2 Two-dimentional function thus could be written as: 
/(mi, ma) = (A-sl+B-bl)-(C-s2+D-b2), where sl(mi) and 
bl(mi) are one-dimentional signal and background functions 
for vector particle, s2(m,2) and 62(7712) _ for pseudoscalar 
particle, A, B, C and D are free fitting parameters. 
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Figure 1: Scatter plots for the processes a) e+e — — > 
(fir) — > i^+i-f - 77, b) e + e~ — » </>r/ — » K + K~ rpr+7r — , c) 
e + e — — > </ir/ — > K + K~ tt + tt~ 7, d) e+e - — > pr; — > tt+tt - 77, 
e) e + e — — > pr/ — » 7T+7T - rpr + 7r — , f) e + e — — » pr/ — + 
Tr + 7r _ 7r + 7r — 7. The dotted lines show the mass ranges used 
for one-dimensional projections. 
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Figure 2: The mass projections for the scatter plot 
m(K + K~ ) vs. 171(77) onto a) K+K~ and b) 77 for the 
reaction e + e — — » </>r/ — > 77. The solid curves show 

the result of the two-dimensional fit, the dotted curves show 
the background contamination. 
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Figure 3: The mass projections for the scatter plot 
m(K + K~ ) vs. m(r/7r+7r - ) onto a) K+K~ and b) 7)7r+7r — 
for the reaction e + e — — + <frq' — > AT+.K" _ r;7r+7r~ . The solid 
curves show the result of the two-dimensional fit, the dotted 
curves show the background contamination. 



backgrounds for the tt + tt , 77, t]it + 'k 



and 



7T 7T 7 



systems are parameterized with first-order poly- 
nomials. The fit ranges are 0.4-1.4 GeV/c 2 , 0.4- 
0.8GeV/c 2 , 0.83-1.12 GeV/c 2 and 0.8-1.2 GeV/c 2 
for 7T + 7r _ , 77, ?77r + 7r~ and 7r + 7r~7, respectively. 
The two-dimensional fitting functions for the scat- 
ter plots are the sum of products of the correspond- 
ing one-dimensional signal and background func- 
tions. The mean values (masses) of the signal func- 
tions are fixed at PDG values [15[ while their widths 
are fixed to the values obtained from the corre- 
sponding inclusive spectra in the data. The sig- 
nificance of the fit is defined as y/— 2\n(L /L max ), 
where Lq and L max are the likelihood values re- 
turned by the fit with signal yield fixed to zero and 
at its best fit value. The signal yields obtained from 
this fit procedure and the significance of the fits for 
all processes are presented in Table [TJ 

The corresponding onc-dimcnsional projections 
together with the results of a two-dimensional fit for 
different reactions are shown in Figs. 2(a,b)-7(a,b). 
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Figure 4: The mass projections for the scatter plot 
m(K+ K~ ) vs. m(7r + 7r - 7) onto a) K + K~ and b) tt + tt~ ■y 
for the reaction e+e — — » </>r/ — + K+K~ 7r+7r~ 7. The solid 
curves show the result of the two-dimensional fit, the dotted 
curves show the background contamination. 
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Figure 5: The mass projections for the scatter plot m(7r+7r _ ) 
vs. m(77) onto a) 7r+7r~ and b) 77 for the reaction e + e _ — > 
pr\ — > 7r+Tr~77. The solid curves show the result of the 
two-dimensional fit, the dotted curves show the background 
contamination. 
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Figure 6: The mass projections for the scatter plot m(7r + 7r _ ) 
vs. m(r?7r+7r — ) onto a) 7r+7r~ and b) rpr+Tr - for the reaction 
e + e — — > pr/ — * TT + TT~rjir + TT~ . The solid curves show the 
result of the two-dimensional fit, the dotted curves show the 
background contamination. 
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Figure 7: The mass projections for the scatter plot m(ir + ir~) 
vs. m(n-+7r~7) onto a) 7r+7r~ and b) 7r+7r~7 for the reac- 
tion e + e~ — > prj' — » tt + tt~tt + tt~ 7. The solid curves show 
the result of the two-dimensional fit, the dotted curves show 
the background contamination. 



From the fitted signal yields, N, we determine the 
corresponding cross sections according to the for- 
mula: 



N 



L By Bp £ 



(1) 



where L is the integrated luminosity (5f6 fb _1 ), 
By, Bp are the branching fractions of the corre- 
sponding decay channels of the vector and pseu- 
doscalar mesons [15[ and e is the corresponding 
detection efficiency. The detection efficiencies are 
determined from Monte Carlo samples where the 
e + e _ — > VP reactions are generated without ISR 
and with an angular dependence corresponding to 
a J p =1" initial state [16j |: 

dN 



d cos 6* d cos 9\ 



sin 2 6 V (1 + cos 2 6* + cos 2<j) V sin 2 67*), (2) 

with the production angle 6* defined as the an- 
gle between the vector meson direction and inci- 
dent e~ beam in the CM frame. The vector me- 
son helicity angle Oy is defined as the polar angle, 
measured in the vector meson rest frame, of the 
positive decay product momentum direction with 
respect to an axis that is aligned with the vec- 
tor meson momentum direction in the CM frame. 
The variable cj)y is the vector meson's positive de- 
cay product azimuthal angle around the direction 
of the vector meson measured with respect to the 
plane formed by the vector meson and the incom- 
ing electron. The generated events were passed 
through a full GEANT [ljj Belle simulation and re- 
construction procedures including the trigger simu- 
lation. The trigger efficiencies, estimated with MC 
samples described above, are about 93% and over 



Table 1: Observed number of events, significance (£) of the Table 2: Total systematic uncertainties for analyzed chan- 

fit, efficiencies and cross sections (cr). nels. 



Process 


N signal 


E 


e, % 


cr, fb 


(£77(77) 


14.6 ±4.3 


8.0 


14.1 


1.1 ±0.3 


(jyq' (rjTT + -K~) 


3.0 ± 1.7 


12.0 


0.917 


2.9 ± 1.6 


<\>r( (n + ir~ 7) 


19.6 ±4.5 


30.0 


5.36 


4.9 ± 1.1 


<f)r]' (comb.) 








4.3 ±0.9 




116.3 ± 20.2 


9.2 


23.2 


2.5 ±0.4 




17.9 ±4.8 


7.9 


3.58 


2.2 ±0.6 


prj' {n + ir~~ 7) 


72.1 ± 15.0 


7.6 


14.3 


3.3 ±0.7 


p?/(comb.) 








2.7 ±0.5 



97% for two-charged-track and four-charged-track 
events, respectively. 

The cross sections before applying radiative cor- 
rections together with the observed numbers of sig- 
nal events, significances of the fit and efficiencies are 
presented in Table [TJ Efficiencies for the processes 
e + e~ — > 0?/(r77r + 7r~) and e + e~ —* pr)'(r)ir + ir~) in- 
clude the branching fraction of r\ — > 77 (39.31 ± 

0.20)% [H. 

The systematic uncertainty on the rj — > 77 detec- 
tion efficiency is dominated by the uncertainty on 
the MC shower simulation in the ECL and other 
material. In order to estimate this uncertainty, we 
compare the ratio of signal yields for the decays 
V ~ > 77 an d V —> 7r 7r°7r° in the data and in Monte 
Carlo simulation. We observe a difference of about 
4%, which is treated as the systematic error in the 
77—5-77 detection efficiency. We assume the uncer- 
tainty in the single photon detection efficiency to be 
2%. The systematic uncertainties due to the exper- 
imenta errors in the branching fractions of the ana- 
lyzed decay channels [15| are 1.3%, 3.4%, 0.5% and 
3.2% for the (prj, <pr]', pr\ and prj' , respectively. The 
systematic uncertainty in the tracking efficiency is 
estimated from D* + — ► D°ir + — ► K~tt + tt + decays 
to be 1% per track. The systematic uncertainty 
from the two-dimensional fit is estimated from vari- 
ations in the number of events with mass values 
and widths floating and fixed. It is estimated to be 
1.5%. The uncertainty in the luminosity measure- 
ment is determined by the accuracy of the Bhabha 
generator, which is 1.4%. The systematic uncer- 
tainty due to trigger efficiency is obtained from 
comparison of the rate of the e + e~ — > <f>i events 
with the one expected from MC simulation. It is 
taken to be 1%. The uncertainty due to limited 



Channel 


Error (%) 


e + e~ - 


-> <fa ~ 


* K + K~n 


5.3 


e + e~ - 


-> 07?' - 


-* K + K~i"fir + TT~ 


7.4 


e + e~ - 


-> (/)!]' - 


-» K + K~ir + Tr~j 


6.2 


e + e~ - 


-> pT] - 


> 7r + 7r _ 77 


5.0 


e + e~ - 


-> prj' - 


-> 7r + 7r _ 777r + 7r _ 


7.0 


e + e~ - 


-> pi] 1 - 


-> 7r + 7r _ 7r + 7r _ 7 


5.9 



Monte Carlo statistics is at most 2% for the pro- 
cess e + e~ — > (jyq' — > K + K~htt + tt~ . The uncer- 
tainty on the charged kaon identification is esti- 
mated by comparing efficiencies of kaon identifica- 
tion in decays D* + — > D°tt + — > K~tt + it + for the 
data and Monte Carlo events. For the uncertainty 
of a charged kaon identification we take the relative 
difference in these efficiencies, which is 0.5% per 
kaon. The systematic uncertainties for all analyzed 
channels are given in Tabled 

To check whether the observed signals are due 
to T(45) decays, we scale the off- resonance signals 
to the on-resonance luminosity, and subtract them 
from the on-resonance signals. The observed num- 
bers of events in the off-resonance data are 1 ± 1, 
1±1, 2 ±1.4, 15.2 ±4.7, 1±1, 7 ±3.6 for the pro- 
cesses 4>rj(n), <fir]' (r]ir + ir~ ) , <f>rj' \it + ir~ 7) , pr](ij), 
pvj '(rjTr + Tt~) and prj' (tt + t:~i), respectively. The re- 
sulting branching fractions for T(4S) — > VP are 
(0.4 ± 0.8) x 10~ 6 , (-0.6 ± 2.8) x 10" 6 , (-0.5 ± 
1.0) x 10~ 6 , (0.8 ± 0.9) x 10~ 6 for the (frrf, prj, 
prj' channels, respectively, which are consistent with 
zero. These results can be expressed as the 90% 
confidence level upper limits [18(, which are equal 
to 1.8 x 10~ 6 , 4.3 x 10~ 6 , 1.3 x 10~ 6 , 2.5 x 10~ 6 
for the (f>r], <frq' , prj, prj' channels, respectively. The 
systematic uncertainties are also taken into account 
for upper limit calculations 

In the light cone approach the authors of 



Refs. [ll|, [l2j gave predictions for the cross sec- 
tions of the reactions analyzed by Belle at s/s = 
10.58 GcV. In Table [3] we present the Belle cross 
sections radiatively corrected according to Ref. [3 
together with theory predictions 



11 



|. The ra- 
diatively corrected cross section can be written as 

(3) 



CT = 



1 



where a is taken from equation (1). The value of 



1 + 5 corresponding to the energy cut of 0.3 GeV is 
equal to 0.809 [l9j | . The BaBar measurement of the 
reaction e + e~ — > <^>'?7 0] is also presented in the ta- 
ble. The corresponding value of 1 + 6 for the BaBar 
energy cut of 0.23 GeV is 0.768. BaBar reports 
that the cross section of the reaction e + e~ — > (prj 
is 2.9 ± 0.5 ± 0.1 fb [6]. The Belle cross section is 
smaller than the BaBar result by about 2.3c. 

From Table [3] we see that in comparison to the- 
ory [ll], [l2| the Belle experimental cross section for 



4>tj is significantly lower, e + e 



pi] 



about 1.8er higher, while e + e~ — > (jnf and e + e~ — > 
prj agree within errors with theory. There is also a 
discrepancy between the data and light cone expec- 
tation in the ratio of the cross section of rj meson 
production together with vector mesons to that of 
rf production. As can be seen from Table [21 the 



light cone approach [ll[ predicts — 



t£n) 



> 1 



r(e+e~ — >PV') 

while this trend is not observed in data. 

The energy dependence of the cross sections 
may have important theoretical implications. In 
Figs.[5Ja-d) we show the Belle data radiatively cor- 
rected according to Ref. [19( together with CLEO 
and BaBar ISR results. The BaBar data were aver- 
aged for y/s values from 2.5 to 3 GeV. We also show 
1/s 3 and 1/s 4 dependences, which pass through the 
CLEO points. In Fig. [SJb) the arrow shows the 
CLEO upper limit and the curves pass through the 
Belle measurement. From Fig. [8] we cannot draw 
any definite conclusion about the energy depen- 
dence of the e + e _ — > VP reactions. 

To summarize, we have measured cross sections 
of the exclusive processes e + e~ — > cjnj, e + e~ — > 
efnf , e + e~ — > prj, e + e~ — > prj where the final 
state includes soft photons with energies below 
0.3 GeV using a 516 fb -1 data sample recorded near 
^fs = 10.58 GeV. The corresponding values of the 
cross sections in femtobarns are: 1.4 ± 0.4 ± 0.1, 
5.3 ± 1.1 ± 0.4, 3.1 ± 0.5 ± 0.1 and 3.3 ± 0.6 ± 0.2. 
The results are compared with theoretical predic- 
tions from the light cone approach 11, . Anal- 



ysis of the energy dependence of the cross sections 
using our results together with those of CLEO at 
y/s = 3.67 GeV [| and BaBar 2.5 GeV for (f>r) and 
2.75 GeV for pr), pi]' Q shows that there is no 
universal energy dependence for these processes. 
The ratios of cross sections of rj meson production 
together with vector mesons to the corresponding 
cross sections for rj meson production are differ- 
ent from the light cone expectation 11 , [13] ■ The 
90% confidence level upper limits on the branch- 




10 , 
s (GeV ) 



10, 
s (GeV ) 



Figure 8: The measured cross sections at yfs « 2.5, 2.75 GeV 
by BaBar, y/s = 3.67 GeV by CLEO and at yfs = 10.58 
GeV by BaBar and Belle for the various processes. BaBar 
measurements are represented by squares, a) e+e — — » {jrrj, 
b) e + e — — > tf>r)'; c) e + e — — » prj; d) e + e — — » pr/' . In the plot 
b) the CLEO upper limit is shown by the arrow. The solid 
lines correspond to a 1/s 4 dependence and the dashed ones 
represent 1/s 3 . 



ing fractions of the T(45) -> VP are 1.8 x 10~ 6 , 
4.3 x 10^ 6 , 1.3 x 10~ 6 , 2.5 x 10" 6 for the <jrr), (pr)' , 
pr], prj channels, respectively. 
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Channel 


cr Belle (fb) 


a [11] (fb) 


a [12] (fb) 


cr BaBar (fb) 


e + e~ — > <prj 


1.4 ±0.4 ±0.1 


3.3-4.3 


2.4-3.4 


2.9 ±0.5 ±0.1 


e + e~ — > cfyq' 


5.3 ± 1.1 ±0.4 


4.4-5.8 


3.5-5.0 




e + e~ — > p?7 


3.1 ±0.5 ±0.1 


2.4-3.1 


2.4-3.5 




e + e~ — > p?/ 


3.3 ±0.6 ±0.2 


1.5-2.1 


1.6-2.3 
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